Observation of a gradient catastrophe generating solitons 
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We investigate the propagation of a dark beam in a defocusing medium in the strong nonlinear 
regime. We observe for the first time a shock fan filled with non-interacting one-dimensional grey 
solitons that emanates from a gradient catastrophe developing around a null of the optical intensity. 
Remarkably this scenario turns out to be very robust, persisting also when the material nonlocal 
response averages the nonlinearity over dimensions much larger than the emerging soliton filaments. 



Introduction In many physical systems propagation 
phenomena are affected primarily by the interplay of dis- 
persive and nonlinear effects. In this context, solitons 
(or solitary waves), i.e. wave-packets that stem from a 
mutual balance between the two effects, account success- 
fully for several phenomena ranging from long-span non- 
spreading propagation and elastic interactions of beams, 
to the coherent behavior of ensembles of particles, e.g. 
ultracold atoms, or coupled oscillators. Studies in this 
field were mainly focused on individual solitons or in- 
teractions between them. However, several solitons can 
emerge at once from breaking of large amplitude smooth 
waves as for instance observed in oceanography 
While theoretical studies indicates the phenomenon to 
be generic [sl, 0], the observation of such multi-soliton 
regime in reproducible lab experiments has been elusive. 

In this letter, we report a lab experiment in optics 
which demonstrates that a fan of non-colliding ID soli- 
tons emerge, owing to a gradient catastrophe (i.e. an 
infinite gradient developing from a smooth input) devel- 
oping around a zero of the field. Specifically we consider 
a dark-like optical beam (i.e., a dark stripe on a bright 
background) and operate, unlike previous experiments on 
dark solitons [I, @] , in a regime where nonlinearity out- 
weighs diffraction (i.e., power of the background largely 
exceeding that needed to trap a fundamental dark soli- 
ton). In this regime, we are able to monitor directly 
the evolution along a thermal defocusing medium. We 
observe the formation of a dark focus point which corre- 
sponds to a gradient catastrophe of the hydrodynamic 
type around a point of vanishing intensity. The infi- 
nite gradient of the hydrodynamic stage is regularized 
by the presence of weak diffraction, which causes the 
appearance of fast oscillations in an expanding region 
(fan), a feature common to the wide class of so-called 
collisionless or dispersive shock waves (DSW) or undu- 
lar bores, investigated theoretically in several contexts 

0, 13, i, m [U, El 0, Ei, 15,16]. 

In our setting, the DSW is essentially composed by ID 
dark soliton filaments, which become manifest after the 
catastrophe point, and maintain fixed parameters (ve- 



locity and darkness) as soon as they emerge [iG^, also 
not exhibiting the rapid decay into vortices characteris- 
tic of shock waves in superfluids [17]. Our scenario turns 
out to be remarkably robust against the nonlocal charac- 
ter of the nonlinearity, and presents also significant dif- 
ferences with DSW resulting from bright disturbances 
[isl . ITqI . |2qI . In the latter case the catastrophe oc- 
curs, indeed, at finite transverse extension, giving rise, 
in 1+lD, to two symmetric fans (connected by a quasi- 
flat background) [H, [2^ , while the relative oscillations 
change dynamically (i.e., dark solitons in the train have 
always slowly- varying parameters upon propagation). 
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FIG. 1: (Color online) Sketch of the experimental setup. 

Experiment Our sample consists of a BK7 glass cell of 
dimensions 1cm x 4cm x 1cm in the X Y and Z (propa- 
gation) direction respectively, containing a solution (con- 
centration C = 0.055mMol/dm^) of Rodhamine-B in 
Methanol. A beam at Aq = 532 nm from a diode-pumped 
CW Nd- Vanadate laser was focused down to a strongly 
elliptical beam (ellipticity 1:30) by means of a cylindrical 
lens of focal length Lf = 100mm and a 20X microscope 
objective. The beam is coupled in the cell at Z = lying 
in the cell mid-plane (sufficiently far from bottom, top, 
and lateral liquid-glass interfaces to avoid the dynamics 
to be significantly affected by boundary conditions). A 
phase mask is placed on the beam path resulting in an 
abrupt change of tt in the optical phase across the line 
X = 0. After the mask we let the beam diffract shortly 
and then we focus it onto the sample, producing an input 
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background bright beam of dimension 600 x 20 /im, onto 
which a dark stripe (with zero intensity in X = and 
hyperbohc-tangent-hke X-profile) is nested. The stripe 
is parahel to the narrower spot-size {Y direction), while 
along X the dark notch of 25/im FWHM sees a quasi- 
constant background due to the 600/im width. We de- 
tect no significant changes along Y over the propagation 
lengths involved, witnessing that our arrangement mim- 
ics a strict 1+lD {X — Z) setting. 

The power coupled into the sample is measured by 
means of a beam splitter in the front of the laser output 
and a silicon detector. As sketched in Fig. [T] a micro- 
scope and a charge coupled device (CCD) camera allows 
us to collect the light scattered in the vertical (Y) direc- 
tion above the cell, so to produce a direct planar image 
of the relevant beam evolution {X — Z plane). In Fig. [2] 
we show the field intensity distribution collected in this 
way at different input powers Pin. 

At low power (4mW) the dark notch diffracts, broad- 
ening in propagation toward positive Z, whereas at 
higher power the nonlinear thermal response of the sam- 
ple counteracts the diffraction leading to exact counter- 
balance (dark soliton formation. Pin — SOmW)^ and sub- 
sequent overall focusing. By further increasing the power 
to Pin = 260m W (Fig. Efc) and Pin = 600m W (Fig. [2]i), 
the dark notch undergoes to a clear focus point (gradi- 
ent catastrophe or breaking point). Beyond such point, 
the beam undergoes non-trivial breaking forming a DSW 
constituted by narrow dark soliton filaments which fill 
progressively a characteristic fan. The number of dark 
filaments in the fan grows with power. This is also clear 
from the measured far-field relative to the output of our 
cell (i.e., after Z = 1 cm of propagation) displayed in 
Fig. [3l Importantly, these data clearly shows that the 
filaments, once seen in the transverse plane, maintain the 
stripe features of the input (similarly to the bright case 
(2q|). not exhibiting any transverse instability or decay 
into vortices characteristic of other superfiuidity and op- 
tical experiments [HI, [22[ , allowing for a description in 
terms of pure I+ID {X — Z) dynamics. 

Theory The dynamics observed experimentally can be 
studied and understood on the basis of the following gen- 
eralized nonlinear Schrodinger (NLS) model 



dz 2 dx'^ 



Sm/j = —i—ei/j^ 



(1) 
(2) 



where the first equation stands for the paraxial nonlin- 
ear wave equation for = A/\/To^ i.e. the beam en- 
velope A = A(Z, X) normalized to peak intensity /q (in 
the experiment /q = Pin/ where A^ is the background 
beam area). The transverse and longitudinal coordinates 
x^z = X/wq^Z/L are scaled to the waist of the input 
dark notch wq^ and the geometric mean L = ^LniL^ be- 
tween the length scales Lni = {ko\n2\Io)~^ and = kw^ 
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FIG. 2: (Color online) Transverse distribution of the inten- 
sity along the cell {x — z plane), as observed from top scat- 
tered light for four different input powers. Superimposed light 
curves (yellow) are retrieved intensity profiles at Z — 0.6 mm, 
while in (d) the right dark (blue) curve is relative to Z = 2.25 
mm. The diffraction fringes aside the central dark notch are 
due to the focusing system before the sample and reflect the 
slight convex wavefront of the bright background. 
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FIG. 3: (Color online) Far-fleld intensity distribution in the 
transverse plane X — Y at the sample output for different 
input powers: the intensity is collected after about 1 meter of 
free- air propagation. 



characteristic of the nonlinear and diffr active terms, re- 
spectively (n2 is the Kerr coefficient that characterizes 
an index change of the local type An = n2|Ap), and 
a = aoL is the normalized attenuation constant. Such 
scaling allows us to highlight the fact that we operate 
in the weakly dispersive case, where the model mimics 
the quantum Schrodinger equation with the smallness 
parameter e = Lni/L = Lnil^d playing the role of 
Planck constant. The normalized refractive index change 
= koLniAn acts as a self- induced potential driven by 
the normalized intensity profile \tjj{x)\'^. The free param- 
eter cr^ measures the diffusion length and gives the degree 
of nonlocality of the nonlinear response. This model de- 
scribes the nonlocal features of the nonlinear response 
with sufficient accuracy regardless of their physical ori- 
gin. Specifically, the model was shown to_give an excel- 
lent description of thermal nonlinearity [2l|, |23| , while it 
allows for a reduction to the integrable semiclassical NLS 
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equation in the local and lossless limit = a = 0. 

The essential physics can be explained indeed by the 
latter limit, for which the outcome of numerical compu- 
tations based on Eqs. (1-2) with input V^o(^) = tanh(x), 
are displayed in Fig. [Hfor two different values of e (pow- 
ers). The initially dominant nonlinearity allows us to 
adopt a description in terms of hydrodynamical vari- 
ables p and u = dxS. This is made by applying the 
WKB transformation ?/^(x, z) = y/ p(x, z) exp [iS{x, z)/ e] 
0, B @, 0, which allows to reduce Eqs. (1-2), at 
lowest order in to the following system written in the 
form of hyperbolic conservation laws 

where q{z,x) = p{z ^ x)u{z ^ x) . Equation (|3]), that rules 
classical ID dynamics of an isentropic gas or shallow wa- 
ter {u being, in this case the velocity of the gas or water, 
and p the gas density or the water level), predicts that the 
dynamics of the input hole in the density p(x,0) = l^/^oP 
produces a gradient of "velocity" whose sign {u turns 
out to be positive for x < 0, and viceversa) is such to 
give rise to compressional waves. Equivalent ly, due to the 
defocusing nature of the medium the central dark region 
has a higher index which draws light inwards. As a result 
the input dark notch experiences a dramatic steepening 
and focusing around its null, which in turn enforces the 
velocity gradient, until eventually a singularity (gradient 
catastrophe) develops at a finite distance, consistently 
with the hyperbolic nature of Eqs. (|3]). Such singularity 
is characterized by crossing of characteristics associated 
with Eqs. (|3]) and become manifest as a vertical front 
in the variable u and a cusp in the intensity p, as dis- 
played in Fig. [Hi (numerical results from Eqs. ([3]) are 
exactly coincident to those shown). However, when such 
a high (virtually infinite) gradients develop, the hydro- 
dynamic description breaks down, and diffraction regu- 
larizes the front through the appearance of an expanding 
region of fast oscillations, i.e. a shock fan. The shock 
fan is progressively filled with non-interacting dark fila- 
ments, whose angle (transverse velocity) increases as the 
relative darkness decreases, which is a universal feature 
of dark (grey) solitons Furthermore the number of 
filaments (solitons) increases for smaller e (higher pow- 
ers). In particular, for 1/e = A^, N integer, one can show 
that the asymptotic state is made by 2A^ — 1 dark soli- 
tons, while no radiation at all is produced [l^. In other 
words, the solitons which are embedded in the input be- 
come manifest at sufficiently large distance, though this 
occurs in a non trivial way through a critical behavior 
characterized by a cooperative initial stage which results 
into the catastrophe focus point. Our experiment pro- 
vides evidence that nonlocality and losses do not quali- 
tatively affect this scenario (at variance with observations 
for bright Gaussian input, which involves a lot of radia- 
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FIG. 4: (Golor online) Numerical simulations of Eqs. (1-2) 
with a = a = and ipo = tanh{x): (a-b) Level plots of IV^I^ 
for s = 0.05 (a) and s = 0.02 (c); (b) snapshot of case (a) 
at z = 4; (d) snapshot of case (c) at the catastrophe point 
z = 0.75: frequency u (solid red curve and intensity (black 
solid curve) compared with the input (blue dashed curve). 



tion instead of solitons [2l|). This is ultimately related to 
the persistence of solitons in the presence of terms that 
break the integr ability. We recall that, in our system, 
the losses and nonlocality are intimately related because 
the index change is due to heating caused by the strong 
absorption of the dye, while the nonlocality arises from 
the intrinsic tendency of heat to diffuse [21]. As shown in 
Fig. O where we report simulations based on Eqs. ([T]|2|) 
with = 1, the overall dynamics is quite similar to the 
local one, except for a slight adiabatic broadening of the 
solitons due to the losses. The post-catastrophe breaking 
still occurs in the form of narrow soliton filaments, which 
model ([IE]) support also in the nonlocal case [2J]. This 
occurs despite the fact that the induced potential Sn that 
trap them becomes, owing to diffusion, a smooth function 
(see Fig. [SJd, red curve) that no longer follows the deep 
oscillations of the intensity as in the local case. Simu- 
lations in transverse 2D with elliptical input reminiscent 
of the experiment confirm the validity of this ID picture, 
allowing us to conjecture that nonlocality stabilizes the 
soliton stripes in the fan against transverse instabilities. 

Importantly, the breaking distance (point in Z of max- 
imum intensity gradient) turns out to be significantly af- 
fected by the attenuation and the degree of nonlocality, 
an issue that we have investigated in detail. In the lo- 
cal and loss- less case, the hydrodynamic limit yields a 
constant normalized breaking distance = 0.75, corre- 
sponding to a physical distance Z5 = z^L that scales with 

— 1/2 

power as P^^ . This is confirmed by numerical solutions 
of Eqs. ([IJ2]) with a = a = 0, performed for s = l/N, 
N integer (i.e. Pin/Pg = A^^, Ps being the fundamental 
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of strong nonlinearity. The post-breaking dynamics gives 
rise to a fan of non- interacting ID soliton filaments. Due 
to weak dispersion such filaments are very narrow, yet 
they are robust against nonlocal averaging over much 
larger widths. 

The research leading to these results has received 
funding from the European Research Council under the 
European Community's Seventh Framework Program 
(FP7/2007-2013)/ERC grant agreement n. 201766. 



FIG. 5: (Color online) Evolution in the nonlocal and lossy 
case: (a) Level plot of intensity (the color scale is adapted to 
compensate for losses); (b) snapshots at z = 7, of intensity 
(black solid curve) and relative index change 5n (red solid 
curve). The input is the dashed blue curve. Here e = 0.05 
(P 0.6W in our experiment), — 1, a — 0.3. The corre- 
sponding physical breaking (catastrophe) distance turns out 
to be Zb — Lzh — 0.8mm. 
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FIG. 6: (Color online) Breaking distance Zb vs. input 
p ower Pi n (log- log scale): (a) local case, numerical results 
(dots) vs. behavior Zh oc P^^^^ (solid line), characteristic 
of the hydrodynamic limit; (b) nonlocal case: experimental 
results (dots) vs. data extrapolated from numerical simula- 
tions (dashed line) performed different values of powers Pin 
(e). Here a = 0.3 was fixed to find the best agreement with 
the data. 



dark soliton power [5[). As shown in Fig. [6^, Z5 ap- 

— 1/2 

proaches the law P^^ for high enough values of TV. A 
similar trend is confirmed by numerical simulations per- 
formed in the nonlocal case by employing the parameters 
of our experiment (we measured by means of Z-scan ap- 
paratus ao = 1.17 mm~^, 77,2 = —2 x lO'^^w? jW), as 
reported in Fig. [6)3. As shown, the expected breaking 
distance agrees reasonably well with the measured data 
except for very low and high powers, the latter show- 
ing a saturation effect not accounted for by our model. 
Here we have used cr as a free parameter, finding the 
best agreement for a 0.3, which is consistent with our 



DtPo Cp\n2\ 
ao\dn/dT\w^ 



1/2 



0.33 [Hj, 



independent estimate a 

based on the values of the parameters for methanol Dt 
IO-W/5, po = 791kg/m^, Cp = 2 x 10^Jkg-^K-\ 
dn/dT = —4 x 10~^i^~^. In summary, we have pre- 
sented the first demonstration of a gradient catastrophe 
occurring around a point of vanishing field in the regime 
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